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We investigate the consequences of an imperfect dark energy fluid on the large scale structure. 

A phenomenological three parameter fluid description is used to study the effect of dark energy on 
the cosmic microwave background radiation (CMBR) and matter power spectrum. In addition to 
the equation of state and the sound speed, we allow a nonzero viscosity parameter for the fluid. 
Then anisotropic stress perturbations are generated in dark energy. In general, we find that this 
possibility is not excluded by the present day cosmological observations. In the simplest case when 
all of the three parameters are constant, we find that the observable effects of the anisotropic 
stress can be closely mimicked by varying the sound speed of perfect dark energy. However, now 
also negative values for the sound speed, as expected for adiabatic fluid model, are tolerable and 
in fact could explain the observed low quadrupole in the CMBR spectrum. We investigate also 
structure formation of imperfect fluid dark energy characterized by an evolving equation of state. 
In particular, we study models unifying dark energy with dark matter, such as the Chaplygin gas 
or the Cardassian expansion, with a shear perturbation included. This can stabilize the growth 
of inhomogeneities in these models, thus somewhat improving their compatibility with large scale 
structure observations. 
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I. INTRODUCTION 

Dark energy is a fundamental component of the nowa¬ 
days standard cosmological model. It would be very 
difficult to explain the set of present days cosmologi¬ 
cal observations without it. Specifically, we refer to the 
luminosity-redshift relationship from observations of su¬ 
pernovae of type la (SNIa) ji [2S, y], the matter power 
spectrum of large scale structure as inferred from galaxy 
redshift surveys like the Sloan Digital Sky Survey (SDSS) 
0 and the 2dF Galaxy Redshift Survey (2dFGRS) |l|, 
and the anisotropies in the Cosmic Microwave Back¬ 
ground Radiation (CMBR) [§. 

Despite of its major importance in explaining the as- 
trophysical data, the nature of dark energy is one of the 
greatest mysteries of modern cosmology. The simplest 
and most popular candidates for it are the cosmologi¬ 
cal constant (see ag. 7]), and minimally coupled scalar 
fields (see e.g. !§I 9L 1 1 tA flTl ). However many other can¬ 
didates were proposed based on hi gh e ner gy p hys ics p he- 
nomenology ( see e.g. [T^, 0, IT~lLT ubL Hit 1 171 iTsl llfll l2fl 
mm ), and many investigations on their possible as- 
trophysical and cosm olog ical signature were undertaken 
^_see_j. g . 

With so many possible candidates it is imperative to 
understand what are the main properties of the dark en¬ 
ergy component that could have specific signatures in the 
astronomical data, and so could help us to discriminate 
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among all these models. 

In a phenomenological approach, dark energy might be 
mainly characterized by its equation of state w , its sound 
speed c s , and its anisotropic stress a (21]. Much effort has 
been put into determining the equation of state of dark 
energy, in an attempt to constrain theories. The equation 
of state determines the decay rate of energy and thus af¬ 
fects both the background expansion and the evolution of 
matter perturbations (see e.g 0]). An equally insight¬ 
ful characteristic of dark energy is its speed of sound. 
This does not affect the background evolution but is fun¬ 
damental in characterizing the behavior of its perturba¬ 
tions. Hence many authors have explored its effect on 
the evolution of fluctuations in the matter distribution 
( see e.g. [Ii| HtJ |u], ^ ) • However, the investigation 
of the effects of the anisotropic stress has been widely 
neglected. 

The main reason for disregarding the anisotropic stress 
in the dark energy fluid might be that conventional dark 
energy candidates, such as the cosmological constant or 
scalar fields, are perfect fluids with a = 0. However, 
since there is no fundamental theoretical model to de¬ 
scribe dark energy, there are no strong reasons to stick 
to such assumption. In fact, dark energy vector field can¬ 
didates have been proposed j43t 144 l45l l4fiL 1471 ] . and these 
have it / 0. Of course, if dark energy is such a vector, 
one might break the isotropy of a Friedmann-Robertson- 
Walker universe. However, as long as it remains sub¬ 
dominant, this violation is likely to be observationally 
irrelevant {48) ]. Once dark energy comes to dominate 
though, one would expect an anisotropic expansion of the 
universe, in conflict with the significant isotropy of the 
CMBR |4!|. But on the other hand there appears to be 
hints of statistical anisotropy in the CMBR fluctuations 
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MISJ.i, 53, ME 3 . 

Recently the possibility of viscous dark energy has 
gained attention J5(| |53, Ha, |53. These models are usu¬ 
ally restricted to the context of bulk viscosity, although 
one could expect the shear viscosity to be dominant [oil . 
One can allow bulk viscosity in a Friedmann-Robertson- 
Walker (FRW) universe, but when the shear is not ne¬ 
glected one has to face the difficulties of an anisotropic 
universe. However, shear viscosity at the perturbative 
level is compatible with the assumption of an isotropic 
FRW background. In fact the anisotropic stress pertur¬ 
bation is crucial to the understanding of evolution of in¬ 
homogeneities in the early, radiation dominated universe. 
Therefore an obviously interesting question is whether 
present observational data could allow for an anisotropic 
stress perturbation in the late universe which is domi¬ 
nated by the mysterious dark energy fluid. 

Motivated by all these possibilities, we investigate if 
the possible existence of an anisotropic stress in the dark 
energy component would result in a specific cosmological 
signature which could be probed using large scale struc¬ 
ture data, and if it would still be compatible with the 
latest CMBR temperature anisotropies and the matter 
power spectrum. 

The article is organized as follows: In section II we 
discuss the parameters describing a general dark energy 
fluid with anisotropic stress. In section III we consider 
dark energy imperfect fluid models parameterized with a 
constant equation of state, sound speed and anisotropic 
stress. We investigate the effects on the late time pertur¬ 
bation evolution, in the integrated Sachs-Wolfe (ISW) ef¬ 
fect of the CMBR anisotropies and on the matter power 
spectrum. In section IV we extend the analysis to mod¬ 
els unifying dark energy with dark matter. We end the 
article with a summary of our findings and conclusions. 


II. DARK ENERGY STRESS 
PARAMETERIZATION 


In its simplest descriptions the dark energy component 
is described fully by its equation of state, defined as 


where p is the energy density and p is the pressure of 
the fluid. If is the four-velocity of the fluid, and the 
projection tensor h^ is defined as h^ = g^ v + 
we can write the energy momentum tensor for a general 
cosmological fluid as 

= pu^Uv -(- ph^v V (2) 

where can include only spatial inhomogeneity. We 
define perfect fluid by the condition = 0. If in ad¬ 
dition the fluid is adiabatic, p = p(p), the evolution of 
its perturbations is described by the adiabatic speed of 


sound c a - This is in turn fully determined by the equa¬ 
tion of state w, 


c 


2 

a 


p W 

— = W -. 

p 3H(l + w) 


( 3 ) 


For an adiabatic fluid, Sp = c 2 Sp. 

In the general case, there may be more degrees of free¬ 
dom and the pressure p might not be a unique function 
of the energy density p. An extensively studied example 
is quintessence !jj. Qjj, jll} . For such a scalar field the 
variables w and c 2 depend on two degrees of freedom: the 
field and its derivative, or equivalently, the kinetic and 
the potential energy of the field. Then the dark energy 
(entropic) sound speed is defined as the ratio of pressure 
and density perturbations in the frame comoving with 
the dark energy fluid, 
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2 
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( 4 ) 


In the adiabatic case, c 2 = c 2 , which holds in any frame, 
but in general the ratio Sp/Sp is gauge dependent. Hence, 
in the case of entropic fluid such as scalar fields, one needs 
both its equation of state and its sound speed as defined 
in Eq. 0 , to have a complete description of dark energy 
and its perturbations. 

However, in order to have an even more general set of 
parameters to fully describe a dark energy fluid and its 
perturbations, besides w and c s , one should also consider 
the possibility of anisotropic stress. This is important 
because it enters directly into the Newtonian metric, as 
opposed to w and c s which only contribute through the 
causal motion of matter j3^ . 

Taking this gene ralization into account, in the syn¬ 
chronous gauge |6Cj, the evolution equations for the dark 
energy density perturbation and velocity potential can 
be written as [6l| 


6 — — (1 + w) | [k 2 + 9H 2 (c 2 — c 2 a )\ ^ + -j 
- 3H(c 2 s -w)S, (5) 


0 = -H{1 - 3c 2 )9 + 5 - k 2 a, (6) 

where h is the trace of the synchronous metric pertur¬ 
bation. Here a is the anisotropic stress of dark energy, 
related to notation of Eq.|0) by (p + p)a = — (kikj — 
4 Sij ) . Basically, while w and c 2 determine respec¬ 

tively the background and perturbative pressure of the 
fluid that is rotationally invariant, a quantifies how much 
the pressure of the fluid varies with direction. 

Generally such a property implies shear viscosity in 
the fluid, and thus its effect is to damp perturbations. 
A covariant form for the viscosity generated in the fluid 
flow is @1 

= C (w/i;a*“ + U t ,;aA“ - U a .J l M „) + ( 7 ) 
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Now the the conservation equations T = 0 reduce to 
the Navier-Stokes equations in the non-relativistic limit. 
Here is the shear viscosity coefficient, and C, represents 
bulk viscosity. Here we set the latter to zero since we 
demand that Ey is traceless. In cosmology we have u^ = 
(1 ,—v } i)/a in the synchronous gauge, and the velocity 
potential 9 is the divergence of the fluid velocity v. It 
is then straightforward to check that the components of 
Eq. ® vanish except in the off-diagonal of the perturbed 
spatial metric. One finds that 

<j=^-( 0-H t ), ( 8 ) 

where Ht is the scalar potential of the tensorial met¬ 
ric perturbations, which in the synchronous gauge equals 
Ht = —/z./2 — 3 ? 7 , where rj is a metric perturbation. From 
the coordinate transformation properties of it follows 
that a must be gauge-invariant, and indeed the linear 
combination 9 — Ht is frame-independent. 

However, the anisotropic stress is not necessarily given 
directly by 9 — H. For neutrinos this term instead acts as 
a source for the anisotropic stress, which is also coupled 
to higher multipoles in the Boltzmann hierarchy. Thus 
the evolution of the stress must, at least in principle, be 
solved from a complicated system of evolving multipoles. 

The approach we will use in this article to specify the 
shear viscosity of the fluid is more in line with the neu¬ 
trino stress than Eq. ©. Following Hu [§3, we describe 
the evolution of the anisotropic stress with the equation 

& + 3H ^=l^ e + -2 + ^ 

Then the shear stress is not determined algebraically 
from fluctuations in the fluid as was the case in Eq. 0 , 
but instead it must be solved from a differential equation. 

This phenomenological set-up is motivated as follows 
JiSl ■ One can guess that the anisotropic stress is sourced 
by shear in the velocity and in the metric fluctua¬ 
tions. Again one must take into account the coordinate 
transformation properties of a , and construct a gauge- 
invariant source term in the differential equation. As 
mentioned, an appropriate linear combination is 9 — Ht- 
Up to the viscosity parameter cA s , this determines the 
right hand side of Eq. ©. In the left hand side there 
appears also a drag term accounting for dissipative ef¬ 
fects. We have adopted a natural choice for the dissipa¬ 
tion time-scale, r ” 1 = 3 H. 

One may then check that Eq. © with w = c^, is = 1/3 
reduces to the evolution equation for the massless neu¬ 
trino quadrupole in the truncation scheme where the 
higher multipoles are neglected 01 (this applies also to 
photons when one ignores their polarization and coupling 
to baryons). In what follows, we will study the conse¬ 
quences of Eq. © for fluids with negative equations of 
state. For w < — 1, one should consider negative values 
of cA s , as was suggested in Ref. (0. So the parameter 
c Sis/(l + w) should remain positive. We will return to 
this in the section III.C. 


Note that the parameterization of Eqs. ©, © and 
© describes cosmological fluids in a very general way. 
The system reduces to cold dark matter equations when 
(w, cl, cl is ) is ( 0 , 0 , 0 ) and relativistic matter corresponds 
to (1/3,1/3,1/3). A scalar field with a canonical kinetic 
term is given by ( w(a ), 1, 0), where —1 < w(a) < 1. With 
an arbitrary kinetic term one can construct k-essence 
models 0 characterized by unrestricted equations of 
state and speeds of sound, (w(a), Cg(a), 0 ), but vanishing 
shear. On the other hand, one should keep in mind that 
the parameterization cannot be completely exhaustive. 
It does not cover, for example, a cosmological fluid with 
anisotropic stress determined by Eq. 0 when c 7 ^ 0. We 
might address the viability of this approximation else¬ 
where, but restrict here to the parameterization given in 
Eq.©. 


III. STRESSED DARK ENERGY FLUID 


We will investigate the effect of dark energy pertur¬ 
bations on the CMBR anisotropies and on the matter 
power spectrum with the simplest assumption that all of 
the three parameters w, cl and cl is are constant. This 
is an accurate description for a wide variety of models 
for which these parameters can be well approximated at 
moderate redshifts by their time-averaged values. 

For the energy content of the universe we use fif, = 
0.044, Qcdm = 0.236, Fide = 0.72 and ho = 0.68. In 
all the numerical calculations we assume a scale invari¬ 
ant initial power spectra and set optical depth to last 
scattering to zero. We normalize the perturbation vari¬ 
ables by setting the primordial comoving curvature per¬ 
turbation to unity. To compare with observations, the 
resulting power spectra must then be multiplied by the 
primordial amplitude of the curvature perturbation. For 
this we employ the same normalization for all the models 
considered in this article. Although we do not search for 
the best-fit models here, we include the WMAP data jQ 
and the SDSS data 0 in to the figures in order to give 
an idea about the viability of the studied models. The 
WMAP error bars include the cosmic variance, which is 
the dominant source of uncertainty for the small Fs. The 
calculations are performed with a modified version of the 
CAMB code |65j . 

In addition to specifying the cosmological parameters 
and the evolution equations ©, © and © , we must 
address the initial conditions in the early universe. For 
the relative entropy between radiation and dark energy 
to vanish, one must impose 


S 

9 


^(1 + w)S r , 

1 

l + 9f^(c2 -c 2 a ) 




■ ( 10 ) 


When there is no inherent entropy in the dark energy 
fluid, i.e. cl = c 2 a , adiabaticity means that the velocity 
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potentials of all cosmic fluids are equal. On the other 
hand, when (? s ^ cl, we see that the above condition 
implies that the dark energy velocity potential is negli¬ 
gible in the early universe, since the relevant scales are 
far outside the horizon, k <C H. However, the condition 
that the relative entropy Eq. m vanishes, when cl = cl, 
would be strictly valid only for the instant that it is im¬ 
posed. Thereby, although is therefore some arbitrariness 
in the choice of initial values, the late evolution is not 
affected by the choice of these values as long as they are 
inside some reasonable region. We use, for all models, 
the initial values 

6 = —(1 + w)6 r , 9 = 9 r , <7 = 0. (11) 

The two first initial conditions are derived assuming that 
c 2 = c 2 and that the relative entropy between dark en¬ 
ergy and radiation together with its first derivatives van¬ 
ish. The third condition says that we set the anisotropic 
stress to zero at very large scales at an early time. Only 
then shear is not generated when c 2 is = 0. Thus fluids 
with vanishing viscous parameter are perfect. Note how¬ 
ever that when the parameter is allowed to evolve with 
time, one can have c 2 is = 0 for an imperfect fluid at some 
stage. 

A. Dark Energy models with — 1 < w < 0 

When w is negative but larger than — 1, the evolution 
of perturbations is determined by the two sound speeds 
of dark energy. The evolution has been analyzed in Ref. 
fsej but without including the anisotropic stress. 

The metric perturbation h is now a source in Eq. ©, 
which tends to draw dark energy into overdensities of 
cold dark matter. However, for large scales the source 
due to velocity perturbations is proportional to — (1 + 
w)(c 2 — w)9/k 2 , and this term can dominate the metric 
source term and drive 6 to smaller values. In fact <5 drops 
below zero when evaluated in the synchronous gauge 1 . 
This happens especially for large sound speeds, since then 
both the friction term in Eq. © and the source term in 
Eq. © are larger (see FIG. ©. Thus 5 gets smaller 
when dark energy begins to dominate. The ISW effect is 
enhanced when one increases the sound speed squared. 

The effect of the anisotropic stress is also to wash out 
overdensities. This is because the metric part of the 
source term in Eq. © turns out negative, and it domi¬ 
nates over the velocity term. Thus a is driven to negative 
values, and in Eq. © it will act to increase the growth 
of 9. This is similar to free-streaming of neutrinos, al¬ 
though for them the effect is relevant at smaller scales. 


1 More accurately, we evaluate the transfer functions of the pertur¬ 
bations. A negative value for the transfer function indicates that 
a perturbation variable acquires the opposite sign to its initial 

value. 



FIG. 1: Late evolution of the dark energy density pertur¬ 
bation and velocity potential for k = 1.3 • 1CP 4 Mpc^ 1 when 
w = —0.8. Solid lines from top to bottom correspond to S, and 
dashed lines from bottom to top correspond to (1 + w)H9/k 2 
when {cl, c 2 vis ) = (0,0), (0.6,0), (0,0.6), (0.6,0.6). The ef¬ 
fect of cl is is to damp density perturbations, which in the 
synchronous gauge is seen as a consequence of enhancing the 
velocity perturbations. 



FIG. 2: Late evolution of the gravitational potentials at large 
scales (k = 1.3 • 10” 4 Mpc -1 ) when w = —0.8 and cl = 0. 
Solid lines are for the case of perfect dark energy and dashed 
for the imperfect case with cl is = 1.0 The upper lines are ip, 
the lower lines are (p. 


Since now cl < 0, the source term ~ H 2 9/k 2 in Eq. © 
inhibits structure growth at large scales. Therefore, as 
the dark energy becomes dominant, the overall density 
structure is smaller when (? vis is larger, and the ISW ef¬ 
fect is amplified. 

It is illuminating to describe the same thing also in 
terms of the Newtonian gauge perturbations. This gauge 
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is defined by the line element 

ds 2 = a 2 (r) [—(1 + 2 (p)dr 2 + (1 — 2ip)dx l dxj\ . (12) 

Here r is the conformal time. We remind that the ISW 
effect stems from the time variation of the metric fluctu¬ 
ations, 


c\ sw 


oc 


dk 

~k 



2 


dr(<j> + ip)j e (kT) , 


(13) 


where tlss is the conformal distance to the last scat¬ 
tering surface and the fth spherical Bessel function. 
The ISW effect occurs because photons can gain energy 
as they travel trough time-varying gravitational wells. 
These wells are in turn caused by matter, since 


— k 2 cp = 47 rGa 2 p 


S +3^(1+ w)6 


(14) 


J|T 


We have indicated with the subscript T that in the left 
hand side variables refer to all matter present, and not 
just dark energy. Note also that the term in square brack¬ 
ets is gauge-invariant. Thus, evaluated in any frame, it 
equals St, the overdensity of energy seen in the comov¬ 
ing frame. During matter domination, St grows in such 
a way that the gravitational potentials stay constant. It 
is then clear that as dark energy begins to take over, the 
gravitational potential \(p\ begins to decay. Contrary to 
expectations from FIG.HJ this decay is not more efficient 
at large scales when there is shear, as shown in FIG. |3 
This is because the dark energy shear influences gravi¬ 
tational wells in such a way that the growth of matter 
perturbations does not slow down as much as in a perfect 
universe. 

However, there is an important twist to the story. This 
is seen in the FIG.|3 where the evolution of the potentials 
4> and ip is plotted at very large scales. At an early time 
the potentials are unequal because of the free stream¬ 
ing of radiation. However, our attention is now on the 
late evolution of the potentials. Due to dark energy, the 
potentials can re-depart from each other at smaller red- 
shifts. This can happen only when c 2 is ^ 0, since 

ip = <f> — 12TrGa 2 (l + w)pcr\T, (15) 


i.e. shear is the difference between the depth of matter- 
induced gravity well and the amount of spatial curvature. 
Since a is gauge-invariant, and we found that it becomes 
negative for dark energy, we can see that shear pertur¬ 
bation drives \tp\ to vanish more efficiently. Thereby we 
find that the effect of shear on Eq. only partly com¬ 
pensates for the effect on ip from Eq. and thus the 
overall ISW from Eq. Cl will be amplified when dark 
energy perturbations tend to smooth as in FIG.^ 

In FIG. m we show the large angular scales of the 
CMBR spectrum when w = —0.8 and the two other pa¬ 
rameters are varied. The upper panel depicts the case 
where the sound speed of dark energy vanishes. Then 




FIG. 3: The CMBR anisotropies for w = —0.8. In the upper 
panel <? s — 0 and in the lower panel c 2 = 1.0. The ISW 
contribution increases with the parameter c^ s : thick lines are 
for c 2 ia = 0, dash-dotted for c 2 ia = 0.001, dashed for c 2 is = 
0.01, dotted for c 2 is = 0.1 and the solid lines for c 2 is = 1.0. 


the pressure perturbation vanishes and the clustering of 
dark energy is inhibited only by the free-streaming effect 
of shear viscosity. Therefore the large scale power of the 
CMBR is increased by increasing c 2 is . In the lower panel 
c 2 = 1. Then dark energy is almost smooth (except at 
the largest scales) even without anisotropic stress, and 
thus we see a smaller effect when c 2 - is increased. 


When c 2 is < 0 the metric and the fluid sources drive 
the perturbations in the same direction, resulting in ex¬ 
plosive growth. Since this would spoil the evolution ex¬ 
cept when c 2 is is tuned to infinitesimal negative values, 
we will not consider such a case here. 
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FIG. 4: Late evolution of the dark energy density perturba¬ 
tion and and the velocity potential for k = 1.3 • ICG 4 MpG 1 
when w = —1.2. Solid lines from bottom to top correspond 
to 5, the dashed lines from top to bottom correspond to 
(1 + w)H8/k 2 when (c 2 , c 2 via ) = (0,0), (0.6,0), (0.-0.6), (0.6,- 
0.6). The effect of c 2 ia is to increase clustering, which in the 
synchronous gauge is seen as a consequence of enhancing the 
velocity perturbations. 


B. Phantom Dark Energy models with w < —1 

When the dark energy equation of state is less than 
— 1, the effect of both the sound speed and of the viscos¬ 
ity parameter are the opposite to the previous case. Now 
the source term in Eq. © has its sign reversed, and be¬ 
cause of that dark energy falls out from the overdensities. 
Similarly, the velocity potential acts now as a source for 
the overdensities. Therefore increasing the sound speed 
will drive dark energy to cluster more efficiently. Now 
the effect of c 2 is > 0 is with the same sign of those of 
the metric sources, and therefore we must consider neg¬ 
ative values for this parameter. Then, if we increase the 
parameter c 2 is /( 1 + w), the dark energy perturbations 
are growing more efficiently, as shown in FIG. This is 
because a is negative, just like in the previous case, and 
again tends to enhance the velocity potential. The cru¬ 
cial difference in the perturbation evolution for imperfect 
dark energy here as compared to the imperfect w > — 1 
case is that more shear in the perturbations will result in 
more clumpy structure in the density of phantom dark 
energy. 

One can again consider the ISW in the terms of the 
Newtonian gauge potentials, Eq. O- The effect is not 
directly seen from the behaviour of S in FIG. 01 partly 
because of the different gauge and partly because the 
anisotropic stress induces a compensation on the other 
gravitional potential and thereby influences also the mat¬ 
ter perturbation. This is shown in FIG. 01 The antici¬ 
pated simple result (that the decay of the gravitational 
potentials is reduced since S is enchanced when there is 


FIG. 5: Late evolution of the gravitational potentials at large 
scales( k = 1.3 • 10~ 4 Mpc -1 ) when w = —1.2 and c 2 = 0. 
Solid lines are for the case of perfect dark energy and, dashed 
for the imperfect case with c 2 is = —1.0. The upper lines are 
<j>, the lower lines are ip. 


more shear) again holds for the sum of the gravitational 
potentials, but considering <p or ip separately reveals the 
intricacy of the fluctuation dynamics due to anisotropic 
stress. Again evolution of </> implies, through Eq. m- 
that the influence of a to dark matter is the opposite 
from dark energy, but on the other hand, evolution of 
the spatial curvature ip implies that the sum <p + ip be¬ 
haves according to the dominating component. Now the 
gravitational well ip grows deeper, because the contribu¬ 
tion from shear in the phantom fluid in Eq. m comes 
with a minus sign. 

In FIG. ©we have plotted the large angular scales of 
the CMBR spectrum when w = —1.2 and the two other 
parameters are varied. The upper panel depicts the case 
that the sound speed of dark energy vanishes. Then 
the ISW effect without anisotropic stress is large since 
dark energy perturbations are nearly washed out. Con¬ 
sequently, the large scale power of CMBR is decreased 
as | c 2 is \ is increased, since the ”anti-viscosity” will then 
amplify perturbations. In the lower panel c 2 = 1. There 
the effect of c 2 already dominates, and we see a smaller 
difference when c 2 is /{l + w ) is increased. 

C. Dark Energy models with c 2 < 0 

For perfect dark energy models without shear, the case 
c 2 < 0 leads to explosive growth of perturbations. This is 
analogous to the behaviour of a simple wave, which has a 
solution ~ e-i(k/c B )t+ik-x, diverging when the sound ve¬ 
locity is imaginary. It is not clear, however, how useful 
the analogy to the sound speed of a simple plane wave is 
to the interpretation of the variable defined by Eq . 0- 
For instance, in the modified gravity context |33j this 
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FIG. 6: The CMBR anisotropies for w = —1.2. In the upper 
panel c a = 0 and in the lower panel c a = 1.0. The ISW 
contribution decreases with the parameter cJi s : thick lines 
are for cj is = 0, dash-dotted lines for cj is = 0.001, dashed 
lines for c 2 is = 0.01, dotted lines for c 2 ia = 0.1 and the solid 
lines for Cy ia = 1.0. 


formal definition does not describe propagation of waves 
in any physical matter. A priori one should not discard 
the possibility c 2 s < 0 without careful deliberation. In 
fact, given a fluid with negative equation of state, one 
would expect, from Eq. 01 , also a negative sound speed 
squared. To get rid of this feature, extra degrees of free¬ 
dom must be assumed to exist in such a way that the 
variable defined by Eq. 0 turns out positive. 

When the generation of shear in the fluid is taken into 
account, the perturbation growth for c 2 < 0 can be sta¬ 
bilized. This is because shear is sourced by the perturba¬ 
tions, and in turn the shear will inhibit clustering. Here 
it is possible to choose the parameters in such a way that 
the dark energy perturbation grows steadily at late times. 
So the ISW effect comes with the opposite sign from the 



FIG. 7: The CMBR spectra for c a = —1.0. The thick line 
is for w = — 1 (unperturbed dark energy). The solid line is 
for w = —0.8 and c 2 ia = 0.5 and the dash-dotted line for 
w = —0.8 and c 2 ia = 1.0 The dotted line is for w = —1.2 and 
c via — 0.5 and the dashed line for w = —1.2 and ct ia = 1.0. 


Sachs-Wolfe effect, which leaves its imprint in the CMBR 
earlier. These effects cancel each other and thus the large 
scale power in the CMBR spectrum is reduced, in accor¬ 
dance with the measured low quadrupole. We show in 
FIG. 0 such a case together with various other choices 
for the other two parameters when c 2 = —1. 


D. A summary 

We summarize the features of dark energy perturba¬ 
tions in different parameter regions in Table ||] A half 
of the parameter space is excluded because divergent be¬ 
haviour occurs, and much of the remaining parameter 
space is degenerated. Even when restricting to simplest 
case where the parameters are kept constant, it seems 
clear that present observational data allows large variety 
of interesting models with non-vanishing shear, c 2 is ^ 0. 
In some parameter regions of Table[I]new features appear 
at observable scales. 

In FIG. M we plot the matter power spectra includ¬ 
ing cold dark matter and baryons for various parameter 
choices. When c 2 , is is varied, the effect occurs only at 
scales much larger than what current observations are 
able to probe. In the spectrum of the total density per¬ 
turbation one would see more pronounced features at 
scales more tantalizingly near the current limits of ob¬ 
servations. However, there is no way to directly measure 
the dark energy density perturbation, and therefore we 
have plotted only the power spectrum of non-relativistic 
matter. 

In FIG. [3 and FIG. 0 it was seen that the shear 
changes the Newtonian gravitational potentials signih- 
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TABLE I: Summary of different parameter regions for dark 
energy fluids. We have indicated with the cases where 
superhorizon perturbations are increased as |c 2 is | is increased, 
and with \ the cases where superhorizon perturbations in 
dark energy are smoothened as c 2 is is increased. We indicate 
by f that the shear perturbation influences significantly also 
the small scale perturbations. 


cantly. Thus one might hope to find a way to study 
whether effects from an anisotropic stress could be mea¬ 
sured by using for example the cross correlation of the 
ISW signal and the large scale structure observations or 
gravitational lensing experiments. However, one should 
keep in mind that we have considered perturbations at 
vast scales. We have found that that fluctuations in an 
imperfect as well as in a perfect fluid with a constant 
equation of state w < —1/3 are confined to super horizon 
scales. This is except for special occasions, in particular 
when the parameter c 2 is negative or when the pertur¬ 
bations behave pathologically due to wrong sign of the 
viscous parameter. For the largest scales, the viscous pa¬ 
rameter determines the evolution of the perturbations. 
It is clear from FIG. 0 and FIG. 0 that the variation 
of the sound speed c 2 has much less effect on the evolu¬ 
tion of perturbations in the limit k —> 0. However, the 
parameter c 2 sets the scale at which the fluctuations in 
dark energy become negligible. For smaller c 2 , there are 
fluctuations at smaller wavelengths. Therefore the shear 
would be best seen when the c 2 is nearly zero or even 
negative. 

The main impact of dark energy anisotropic stress on 
observations seems to be the modification of the CMBR 
at very large scales, from which it would be very dif¬ 
ficult to unambiguously detect. However, this changes 
when one considers the perhaps better physically moti¬ 
vated situation where the parameters w, c 2 and c 2 is are 
allowed to evolve in time. 


IV. IMPERFECT UNIFIED MODELS 

The Chaplygin gas @ is a prototype of a unified 
model of dark matter and dark energy [6§| . In such 
models a single energy component accounts for both the 
dark matter and dark energy. Thus this component 
must resemble cold dark matter in the earlier universe, 
whereas it should exhibit large negative pressure nowa¬ 
days. These models are, however, problematic because of 
the suppression of structure formation by the adiabatic 
pressure perturbations E1IIIE3. A solution for this 
problem has been based on the observation that due to 



FIG. 8: The total matter power spectra when c 2 is = 1.0. 
The thick line is for w = —0.8 and c 2 = 1.0. The dash-dotted 
line for w = —0.8 and c 2 = 0, the dotted line for w = —0.8 
and c 2 = 1. The dashed line is for w = —1.2 and c 2 = 0, the 
solid line for w = —1.2 and c 2 = 1. 


entr opy , the sound speed is not necessarily the adiabatic 
one [HEtzlIzl In the so called silent quartessence model 
entropy perturbations cancel the effect of the adiabatic 
sound speed ^73j. 

The modified polytropic Cardassian expansion 0 
(MPC) in the fluid interpretation (7^ provides a gen¬ 
eral parameterization which encompasses a wide variety 
of unified models. For the MPC case, one can write the 
energy density as a function of the scale factor as 

p = [Aa 3q( ' l '- 1) + Ba~ 3<1 }«. (16) 

The exponents q, v are given as parameters, and the con¬ 
stants A, B have the appropriate mass dimension. This 
parameterization is equivalent to the New Generalized 
Chaplygin gas (pgj. When v = 2, one gets the Gener¬ 
alized Chaplygin gas ,77] where q can vary, and setting 
further q = 1, one is left with the original Chaplygin gas 
§37]. On the other hand, when the parameter v can get 
any values and q is set to q = 2, the Variable Chaply¬ 
gin gas Q is recovered. Finally, when q = 1 and v is 
arbitrary, one has the simplest version of the Cardassian 
expansion which reproduces the background expansion 
of a universe with standard CDM and dark energy with 
w = v- 2 [7|. 

We will consider here effects of shear in models where 
the unified dark density is defined by Eq. (IlGII . In 
the above references, various theoretical origins for such 
density ansatzes have been proposed, but we will not 
consider here whether these necessitate incorporation of 
shear in the linear order in cosmology. Previously Car¬ 
dassian expansion has been studied in the modified grav¬ 
ity context [33^ . and there it was shown that an effec¬ 
tive anisotropic stress can appear in the late universe. 
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With certain assumptions for the modified gravity, the 
cold dark matter density perturbation generates a shear 
perturbation algebraically determined from the density 
and velocity fields of the matter interestingly similarly 
to the case motivated by the covariant generalization 
of the Navier-Stokes equation 0 . Then the resulting 
matter power spectrum is in better accordance with ob¬ 
servations than in the standard adiabatic case, but as 
the anisotropic stress affects the gravitational potentials 
and thus enhances the ISW effect, the CMBR spec¬ 
trum will then restrict the allowed parameter space strin¬ 
gently. However, we will here consider the anisotropically 
stressed fluid as parameterized by Eq. 0 . As expected, 
our results will be different from the modified gravity 
approach of Ref. ■ 

The equation of state for the unified fluid described by 
Eq. lUlifl is 


t/Aa 3q ^~^ 

W = - 

Aa 3 ®!"-!) +Ba~ 3q ’ 

and it follows that 


c 


2 

a 


w [1 — vq + w(\ — q)\ 

1 + W 


(17) 


(18) 


When both q and v are equal to one, the model is equiv¬ 
alent to cold dark matter and a cosmological constant. 
When either q or v is smaller than one, c 2 will be nega¬ 
tive in the late universe, and when either q or v is greater 
than one, the c 2 will stay positive and grow in the late 
universe. For instance, if v = 1, the asymptotic value is 
c* = 1 - <}■ 

In the adiabatic case then c s = c 0 , but for the 
silent quartessence one imposes a special condition on c 2 , 
namely that the pressure perturbation vanishes in the 
synchronous gauge 0 , Z3j • Here we consider the case 
that all the three sound speeds are equal in magnitude, 
including the viscosity parameter as it appears in Eq. 0 . 
Thus (? vis = |c 2 | and c 2 = ( :2 a . This seems a natural gen¬ 
eralization of the characteristics of better known cosmic 
fluids, i.e. neutrinos. Then, in the terminology employed 
here, the fluid is adiabatic but imperfect 2 . For compari¬ 
son, we include also results for the silent quartessence, as 
an example of an entropic but perfect fluid. In addition, 
results for the adiabatic and perfect model are shown. 

We plot the results for the CMBR spectrum in FIG. 
El and for the matter power spectra in FIG. 03 and FIG. 
ni! Here the matter power spectra include all the compo¬ 
nents in the energy budget, since one cannot distinguish 
dark matter from the unified model. In all the figures, 
the adiabatic case is shown with dash-dotted lines, the 
entropic case with dashed lines and the imperfect case 
with solid lines. 


2 The presence of anisotropic str ess does not lead to generation of 

entropy at the linear order t80l . 




FIG. 9: The CMBR anisotropies for the MPC model with 
q = 1.0. The dasli-dotted lines are for the adiabatic case, the 
solid lines for the same model with the shear included, and 
the dashed lines correspond to the silent case. In the upper 
panel v = 1.1. In the lower panel v = 0.9 for the silent and 
the imperfect case. The adiabatic model has v = 0.994, which 
already gives disproportionately large ISW effect. 


The effect of shear is to stabilize the perturbations. 
When q or v is greater than one, the adiabatic pres¬ 
sure tends to drive the density perturbation to oscillate. 
However, including the anisotropic stress will remove the 
oscillations, since the damping effect of shear compen¬ 
sates steep gradients. The overall suppression of growth 
is alleviated, but not as much as in the silent model. 
When either q or v is smaller than one, the adiabatic 
pressure would drive the density perturbation to a very 
fast growth. However, as discussed in the previous sec¬ 
tion, the shear viscosity eliminates this driving force. 


V. CONCLUSIONS 

In this article we have investigated the effects of an 
anisotropic stress in the dark energy component on large 
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FIG. 10: The total matter power spectra for MPC model 
when v = 1.0 In the upper panel v = 1.01, and in the lower 
panel v = 0.999. 


scale structures. 

We have parameterized the dark energy component 
with three variables. The equation of state determines 
the decay rate of dark energy, and the sound speed char¬ 
acterizes the evolution of its fluctuations. These two 
were treated as independent parameters, thus account¬ 
ing for possible entropy in the fluid. In addition we al¬ 
lowed for shear viscosity in the linear order. We discussed 
the possibility to apply a Navier-Stokes type viscosity to 
determine the additional degree of freedom for dark en¬ 
ergy fluctuations, the amount of shear viscosity, but we 
adopted the parameterization utilizing a viscosity param¬ 
eter (? via , motivated by the fact that it seems to generalize 
the familiar and well understood cosmological fluids in a 
natural way @]. 

Using this phenomenological three parameter fluid de¬ 
scription we investigated the effect of an imperfect dark 
energy fluid and of unified dark matter and dark energy 




FIG. 11: The total matter power spectra for MPC model 
when v — 1.0. In the upper panel q = 1.01, and in the lower 
panel q = 0.999. 


models on the matter power spectrum and on the CMBR 
temperature anisotropies. For most models we find that 
free streaming effects tend to smooth density fluctua¬ 
tions. However, there are some exceptions, described 
below. 

In dark energy models where — 1 < w < 0, we found 
that increasing the anisotropic stress results in a swifter 
decay of dark energy overdensities, which is seen in the 
CMBR spectrum as an amplification of the ISW effect. 
The opposite occurs in the case of phantom dark energy 
(w < —1), for which the anisotropic stress supports the 
growth of overdensities and thus reduces the ISW effect. 
However, the impact of anisotropic stress on the CMBR 
spectrum can be closely mimicked by varying the sound 
speed of dark energy. This makes it difficult to distin¬ 
guish between these two fluid properties. 

In addition, we found that negative sound speeds are 
also consistent with observations, if shear viscosity is in- 
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eluded. The situation that the pressure perturbation 
(evaluated in the comoving frame) is of the opposite sign 
than the density perturbation, is formally unproblemat- 
ical to define, but when c^ is = 0 it will exhibit unlimited 
growth of density fluctuations. However, when (? vis > 0 
this does not occur. For a suitable choice of parameters 
a low amplitude for the CMBR quadrupole is produced, 
in accordance with observations. 

In models unifying dark matter and dark energy ex¬ 
tended with shear, it is found that the anisotropic stress 
can stabilize the effect of the adiabatic pressure pertur¬ 
bation, thus slightly improving the compatibility of these 
models with large scale structure observations. It re¬ 
mains to be seen how one can loosen the constraints by 
allowing for an anisotropic stress. Our main objective 
here was to use these models as examples of dark energy 
with evolving w, c s and c v i s . The conclusion taken is 
that, in contrast to the simplest fluid models with con¬ 
stant w, c s and c V i S , in specific scenarios the shear stress 
can have consequences distinguishable with present ob¬ 


servational data. 

In general, we found that anisotropic perturbations in 
dark energy is an interesting possibility which is not ex¬ 
cluded by the present day observational data. Further¬ 
more, we found that the CMBR large scale temperature 
fluctuations, due to the the ISW effect, are a promising 
tool to constrain the possible imperfectness of the dark 
energy component. Even when the anisotropic stress can¬ 
not be directly measured, it can still bias measurements 
of other parameters, for instance the dark energy speed 
of sound or its equation of state. 
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